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We propose an explicit formula for an entanglement measure of pure multipartite quantum states, 
then study a general pure tripartite state in detail, and at end we give some simple but illustrative 
examples on four-qubits and m-qubits states. 
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I. INTRODUCTION 

One of the unsolved problems of modern quantum the- 
ory is the quantification of multipartite state entangle- 
ment 0, 0] ■ This is a task that is directly linked to math- 
ematics, such as linear algebra, geometry and functional 
analysis. The definition of separability and entanglement 
of a multipartite state was introduced in |^ Q follow- 
ing the definition for bipartite states, given in 1989 by 
Werner Jsj. Eventually, quantitative measures, such as 
the entanglement of formation and concurrence were for- 
mulated for bipartite systems |^. In recent years, there 
have been attempts to find an entanglement measure 
for qubit-qudit states B, B. lid and for multipartite 
states, i.e., in [H El ETIM l3 El- To exemphfy, 
upper and lower bounds for the quantum relative en- 
tropy of entanglement of a multipartite systems in terms 
of the bipartite entanglements of formation, distillation, 
and quantum entropy of various subsystems are derived 
in Measures based on the Schmidt rank are pro- 

posed in E3] and on local unitary and filtering processes 
in Furthermore, in 0|, a very useful tool to detect 
entanglement, called entanglement witness, is generalized 
to multipartite states. The tool is a consequence of the 
Hahn-Banach theorem which states that for any convex, 
compact, vector set S, ii p ^ S, there exists a hyper- 
plane that separates p from S. However, to find such an 
operator, even in case of tripartite state, is a formidable 
task. None-the-less, quite impressively, Acin et al. man- 
aged to construct a witness operator for a class of mixed 
tripartite states Ell- 
in a recent paper, Partovi proposes an entanglement 
measure based on generalized Schmidt-decompositions of 
a state E3- I'^ essence, his measure gives the difference 
between the minimum entropy of the separable state hav- 
ing the same statistical marginal distribution as the state 
being characterized, and the entropy of the state itself. 
To find the minimal difference, Partovi employs a suc- 
cessive Schmidt-decomposition of the state. The mea- 
sure quantifies the logarithmic "quantumness" of a state 
irrespective of the type of entanglement is possesses. 



In this paper, we propose another measure of entan- 
glement for arbitrary, pure multipartite states. Inspired 
by the work in j20j. we give an explicit expression for 
such a functional. Our method is based on the joint 
relative-phase properties of a multipartite quantum sys- 
tem Q = Qi ® Q2 ® ■ ■ ■ ® Qm on a Hilbert space 
Hq — <8i Hq^ (81 • • • (8) Hg„ expressed by a posi- 
tive operator value measure (POVM) Ag on Hq. The 
POVM is constructed by taking the m-fold tensor prod- 
uct of the subsystems' corresponding POVMs. We have 
already discussed, in detail, our degree of entanglement 
for a bipartite state in |2lL |23| , so here we will only give 
examples for multipartite states. 



II. ENTANGLEMENT FROM A 
RELATIVE-PHASE POVM 

A general and symmetric POVM in a single iV„ 
dimensional Hilbert space 7iQ„ is given by 



where \ku) are the basis vectors in 7ig„ and 
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The POVM is a function of the iY„(iV„ - l)/2 relative 
phases (<y5i„,2„, • ■ • , <<3i„,Ar„, (<32„,3„, ■ • ■ , 'Pn^-i,nJ- 

It is now possible to form a POVM of a multipartite 
system by simply forming the tensor product 

^Si (VQi;fei,ii) «) • • ■ ® ^Q^i^Q^;k„,,l„J, (3) 

where, e.g., fQuki.ji is the set of POVM relative 
phases associated with subsystems Qi, for all ki,li = 
1, 2, . . . , A^i , where we need only to consider when li > ki 
due to We can now recast this POVM, expressed 
in local properties, in terms of the relative-phase sums 
4>ki,h,...,k^,i^ = J2'u=i Vk^,i^- Note that if, e.g., U = ky, 
then the term (pk^,k^ vanishes from the sum due to 
From Aq we can form an associated real function ex- 
pressed in H^^iNuiNu ~l)/2 ~ M linearly independent 
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relative-phase sums: 

Tr ^/5A(0^(i) ^^(i)^ ^^(1)^^(1), . . . 1 4>^.iJ^n liM) j^ijvi) ,(4) 

where p is the state density operator acting on the com- 
posite Hilbert space Hq. Next, we define to what ex- 
tent the density operator depends on the particular joint 
relative-phase sum ;(i) , e.g., 

^k['\i['\...MJ.\i'^' " 

where V must be expressed in the relative-phase sum pa- 
rameter 0^(1) ^(1) ^(1) ^(1), but the particular choice of 

the remaining M — 1 linearly independent relative-phase 
sum parameters is inconsequential for the absolute value 
of the integral. The coefficients 7fc„,i„,...,fc„,;„ (where, 
here, and in the following, we will omit the superscript on 
the indices) are proportional to the Fourier components 
of the joint relative-phase distribution. Now, let us intro- 
duce the following index operator to connect the notation 
using the subsystem indices, and the notation using the 



joint-system index running from 1 to N1N2 ■ ■ ■ N,n- 
n(fci,;i,A;2,^2, • • • ,k„i,lm) = (ki - 1)N2---N,n 

+ {k2 - 1)N3 ■■■N,n + --- + ik„^-l - l)N,n + k^, 

{h - l)N2 - ■ ■ N,r, + {I2 - l)N:, - ■ ■ Nra 

+ --- + {lrn-l-l)N^+l^. (6) 

Note that the index operator generates two indices 
based on the set {ky} and {lu}^ respectively. Evalu- 
ating the Fourier components, one finds, not surpris- 
ingly, that 7fei,ii,...,fc„,u = 27r|pn(fci,;i,...,fc„,i„)l- That is, 
to each relative-phase sum there is an associated joint- 
system density matrix coefficient. We now define an in- 
dex permutation operator operating on any function 
f{ki,li, . . ..kmjni) by 

^li ■ ■ ■ I kj, Ij, . . . , km, Im) = 
f {k\ , ^1 , . . . , /Cj , /j , . . . , km : ^m) 
~f{k\, ll, . . . , Ij, kj, . . . , km, Im)- (7) 

Using this operator we can generalize our earlier results 
for bipartite systems 0, |23|. We form an entangle- 
ment function by summing the absolute difference be- 
tween pairwise relative-phase sums. The function is given 

by 



I 

ATi TVi N2 N2 N3 Af„ 

^(p) = Yl E •■• \P2\pmiMM,...,k..-i.i^-uk^.i^)\\^ + ■ ■ ■ 

li>ki ki — 1 l2>k2 k2 — l k^—l^ — l km—lm — ^ 

+N2 Y Y X! X! X! ■" X! |Pm|Pn(fci,;i,fc2,i2,...,fc™-i,i™-i,fc™,U)lP 

/m_l>fem_l fcm_i = l im>fe„,A/'2 fcm = l ki=lt=l fc„_2=im-2 = l 

Afi Ni N2 N2 N3 N3 Ni 

J2 J2 J2 11 11 11 H ■■■ H {|P2|P3|pn(fci,h,fe2.i2,fe3.i3,...,fe™-i.i™-i.fc„,;„)ll'|} + --- 

/l>fcl ki — l l2>k2 k2 — l l3>k3 /C3 — 1 «4— J4 — 1 kjn—lm — ^ 

N,„-2 JV,„_2 N„^ JVi Wm-S 

+•^3 H I] £ £ I] H {\Prn-l\Pm\pn{kuh,k2M,...,k^-l,l,r,-lM^)\f\} 

',7i-2>fcm-2 fcm-2 = l I m> km km = i ki =li = l im_3=j„_3 = l 

+ ...+ 

+^r,r-l Y 11 ■ Y H {|P2|P3|---|Pm-l|Pn(fci.ii,fc2,i2...,fe^-2.U-2,fe™-l,U_l,fe^,U)ll'---|ll} 

h>ki fci=l i„i-l>fe„i-l fcn,-l = l km=lm = l 

N2 N2 Nm JV,„ Ni 

+ ...+Nm-lY H H H H {|P3|P4|---|^m|pn(fei,;i,fc2,i2,/c3,i3,...,fc™-2,;™-2,fe™-l,i„.-l,fc„,;„)lP---|||} 

i2>'C2 fe2 = l /m>/c„, fc„j = l fci=/i = l 

Ni Ni Nm Nm, 

+AfmY E--- E E{|P2|P3|---|P™|Pn(fe„;„.2.Z2.....fc_.A.-.,fc.,U)ll'---|ll})^- (8) 
/l>fclfcl = l i„>fe„ fc„=l 

I 

This is our central equation. It looks messy, but has a rather logical inner structure. The factors 
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are normalization factors, and they should not be con- 
fused with the space dimensions N^. The first sums, 
where two of them are written explicitly (with normal- 
ization factors M2) on the right hand side of |(SJ), only 
contributes the respective subsystem's bipartite entan- 
m 



glement. There are 



TO(m — l)/2 ways to se- 



lect two systems out of m without respect to ordering. 
The two terms explicitly written above sums the bipartite 
entanglement contribution between systems Qi,Q2 and 
Qm-i, Qm, respectively. For the systems {Q„, Qy}, there 
are Nu{Nu — l)Ny{Ny — l)/A ways to select one each of the 
relative phases of system Q„ and Qy. Because the other 
system's coefficients can be chosen arbitrarily among the 
diagonals, there are HJLiNj / (NuNy) number of relative- 
phase sums and differences involving fc^ , ^„ , ky and ly . 
Our permutation operator subtracts the relative-phase 
difference from the relative-phase sum, so by including all 
bipartite combinations, the bipartite entanglement of the 
joint system is taken care of. Next, we add the tripartite 
entanglement (contained in the sums with normalization 

m 



factors As). There are 



tripartite combinations. 



and for every choice {Qu, Qv, Qiu}, where u < v < w, 
there are Nu{Ny -l)Ny{Ny- 1)7V^(7V^ - l)/8 combina- 
tion of system relative phases. For each combination, we 
can sum all three relative phases, sum the first two and 
subtract the third, etc. To form differences of all combi- 
nations, we use both the permutation operators P„ and 
Pyj. Hence, we get 4 = 2^~^ contributions within the first 
curly bracket in ((SJ, above. For each choice, the other 
systems indices can be chosen in ITJL^Nj / (NyNyNyj) dif- 
ferent ways. For the quadripartite contribution we pro- 
ceed in the same way. For every choice {Qu,Qv, Qw, Qz}, 
where u < v < w < we use the permutation operators 
Py, Pyj, and Pz- We get 8 = 2'*"^ contributions inside 
the corresponding curly brackets, each being a sum of 
U'p^iNj / (NuNyNyjNz) terms. The sum proceed in this 
fashion until the m-partite entanglement contributions 

are to be added. There is only one way (( I =1) 

to chose all subsystems, and we insert m as index in our 
permutation operator, we use the permutation operators 
P2, P3, . . . , Pm- (We do not permute ki and h.) In all, 
we get 2™~^ terms inside the curly brackets of the last 
sum in ©, above. These terms represent all the possi- 
ble relative-phase sums and differences between all the 
m-systems, so there are no further terms. 

From our definitions, it is clear that for any product 
state 



where pfc„,i„ is the indicated density matrix coefficient of 
system u. In this case, one gets Pu|pn(fci,ii,...,/c„,;„)l = 
for any u and any set of indices fci, Zi, . . . , km, Im- Hence, 
our entanglement function T(p) — for any tensor prod- 
uct of m density operators. For entangled states, the 
function is not invariant to local unitary transformations. 
In analogy with our definitions for bipartite states, we de- 
fine our measure of entanglement Tsup, where sup refers 
to the supremum of T under all possible local unitary 
transformations. 

Let us now write out and use lO in a few explicit cases. 
The degree of entanglement for a Hq^ ® Hq^ bipartite 
states is given by 



Ni Ni N2 N2 



m - (-^2 E E E E ii/'c^.-i 



)JV2+fc2,(/l-l)Af2-|-/2l 



ll>ki ki = l /2>fc2 ^2 = 1 



|P(fcl-l)Af2 + i2,(il-l)A'2 + fc2 IIQ1Q2 



(10) 



This special case has already been discussed in detail in 
21, 22], and we have shown that the equation coincides 
with the concurrence 6] for pure bipartite states in 2 (g) 2 
(provided that one sets JV2 — 2) and with generalized 
concurrence measures in 2 (g 3 dimensions |3j B 13 ■ 

Note that our measure sums all the state's entangle- 
ment. That is, although, e.g., a state's bipartite entan- 
glement between subsystems Qi and Q2 cannot be used 
simultaneously neither with its bipartite entanglement 
between subsystems Qi and Q3, nor, e.g., with its tri- 
partite entanglement between subsystems Qi, Q2, and 
Q3, all contributions are added in our measure. That is, 
our measure characterizes the entanglement contained in 
a state, but in general the measure exceeds the usable en- 
tanglement. However, by looking at the various terms in 
the sum, the usable entanglement can be extracted as the 
measure is composed of sub-sums containing the bipar- 
tite Qi and Q2 entanglement, the bipartite Qi and Q3 
entanglement, the tripartite Qi, Q2, and Q3 entangle- 
ment, etc., as can explicitly be seen in below. Also 
note that our measure sums the possible cooperative en- 
tanglement. That is, if some subsystems are ignored, or 
the information contained in a subsystem is lost, then 
the ensuing state's entanglement is in general lower than 
what our measure predicts. We shall give a concrete ex- 
ample of this in Sec. IIVL below. 



Pn{kuh,k2,l2,---,k^,lm) — PkiJiPk2,l2 ■ ■ ■ Pk,-n,l,-, 



(9) 



III. TRIPARTITE ENTANGLEMENT 

The degree of entanglement for a H-q^ g) 'Hq2 € 
tripartite state is given by 
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Afi Afi N2 N2 N3 

m - mT. E E E E \\Pn{ki,h,k2,hM,h)\ ~ \Pn{ki,h,hM,k3,h)\\QiQ2 

h>kl ki = l l2>k2 k2=l k3=l3 = l 

JVi Ni N3 N3 N2 

+ E E E E E \\P^{ki,h,k2,l2,k3,l3)\~\P^{ki,li,k2,l2,l3M)\\%.lQ3 

h>ki ki = l l3>k3 k3 = l k2=l2 = i 
N2 N2 N3 N3 Ni 

+ E E E E E \\Pn{kuh,k2d2,k3,h)\ ^ \Pn{ki,hM,h,h,k3)\\'Q2Q3^ 

l2>k2 k2 = l l3>k3 k3 = l ki=li = l 

Ni Ni N2 N2 N3 N3 

+■^3 E E E E E {\\\Pniki,h,k2,l2,k3,h)\ - \Pnikuh,k2.h,h,k3)\\'gAQ2Q3 

h>ki ki = l l2>k2 k2 = l l3>k3 k3 = l 




Let us now give two concrete examples of this mea- 
sure for some three-qubit states. In the three-qubit space 
there exist two classes of states, inequivalent under local 
operations and classical communication (LOCC), called 
I^'ghz) and j^-w) states. They are, e.g., I^'ghz) = 
(|000) + |lll))/\/2 and l^-w) = (|001> + |010> + |100>)/V3. 
For these states, we have 




and 

r(pw) - (-^2(|P3,5|q,q, + |P2,5|g,Q3 + \P2.3\%^Q,)) ' 

^ V 3 ' 

Here, we see that the normalization factors must be re- 
tained (or, possibly be chosen with particular relative 
weights) in order for the entanglement measure to make 
sense for general states. Because the GHZ- and the W- 
states belong to different equivalence classes, their rela- 
tive entanglement weights are not obvious. This issue is 
tied to the, still open, question about minimum reversible 
entanglement generating sets 0, 0, 0, 0, . 

IV. BEYOND TREE-PARTITE QUBIT 
ENTANGLEMENT 

Next, let us look at an interesting four-qubit state 
l^-i) = (|0, 0, 0, 1) + |0, 1, 0, 0) + |1, 0, 1, 0) + |1, 1, 1, l))/2. 
Our measure of entanglement of this state is 

rc^i) = i^^2{\P2^5\Q2Q, + \pns6\Q2Q,) 

+A^3(|P2,16|qiQ2Q3 + |P5,1i|qiS2Q3 
+ \P2m\q,Q3Q^ + \P5m\q,Q3qJ)^ 



(11) 

I 

The state has both bipartite and tripartite entanglement. 
In order to use the bipartite entanglement the parties 
possessing the different qubit subsystems must cooper- 
ate. If, e.g., qubit 1 and 3 are measured in the standard 
basis, the result is either two zeros or two ones. If this re- 
sult is communicated to the parties holding qubit 2 and 
4, (that is, we perform a LOCC, optimal for bringing 
out the bipartite entanglement) the remaining two-qubit 
state will be in (a known) pure EPR-state. If, on the 
other hand, if we simply trace out subsystems Qi and 
Q3, (or measure qubit 1 and 3 but keep the result se- 
cret) , then the remaining state is in an an equal mixture 
of the EPR-states, and this state is directly separable. 
This means that if the different parties do not cooperate, 
the state's bipartite entanglement in subspace 1 and 3 
vanishes. 

The entanglement of the state |^'2) = (|0, 1,1,0) + 
|1, 0,0,1) -I- |0, 1,1,1) + |l,0,0,0))/2, on the other hand, 
is given by 



r(f 2) = (AA3|p7,9|Ls.S3)^ = 




That is, the state has only tripartite entanglement and no 
bipartite entanglement. To arrive at the result, we note 
that a unitary transformation t/4 local to Q4 can trans- 
form the state into, e.g., J74|^'2) = (|0, 1,1) + |1,0,0) (g) 
|0))/V2 for which one finds the supremum of F. In this 
case, the state's entanglement is the same whether or not 
the person in possession of qubit 4 cooperates or not. 

As a last example, consider a a m-qubit density oper- 
ator given by a mixture of the two orthogonal, so called, 
m-Cat states 

l^-Cat) = 4s(|0l'02,...,0m) + |ll,l2,...,lm)) 
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. Then, our degree of entanglement gives 

r(fcat) = (-^™IP1,2™Iq,s,...s„)' 

V. CONCLUSIONS 

In conclusion, we have proposed an entanglement mea- 
sure for pure multipartite quantum states. The measure 
directly detects product states (it is zero for such states), 
and quantifies the entanglement of any pure state up to 
the bipartite, tripartite, . . . , m-partite normalization co- 
efficients. Since it is not possible to use, nor convert, the 
entanglement in states with incompatible entanglement 
classes such as GHZ- and W-states into each other, it 
may not be meaningful to specify the coefhcients rela- 
tive to each other. Rather, from an operational point of 
view, it seems more meaningful to specify each type of 
entanglement separately, e.g. in a system composed of 
four subsystems Qi, Q2, Qs, and Q4, it is meaningful to 
discuss, separately, the bipartite entanglement between. 



e.g., systems Qi and Q2, and Qi and Q4. We do not 
see how the bipartite entanglement, in an operational 
sense, could (or should) be compared to, e.g., the tripar- 
tite entanglement between Qi, Q2, and Q4. Our measure 
sums all contributions to quantify the state's entire en- 
tanglement, but, as just indicated, from an operational 
viewpoint, it is rather the the sum's various contribu- 
tions that have a well defined operational meaning. This 
is in contrast to, e.g., Partovi's measure that is a 
minimum entropic distance measure between the state 
and a separable state with the same statistical marginal 
distributions. In this sense, Hossein Partovi's measure 
is a better measure of the "quantumness" of the state, 
while our measure emphasize the state's usefulness as a 
quantum information carrier. 
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